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ABSTRACT 



We use Hirota's method formulated as a recursive scheme to construct complete set of 
soliton solutions for the affine Toda field theory based on an arbitrary Lie algebra. Our 
solutions include a new class of solitons connected with two different type of degeneracies 
encountered in the Hirota's perturbation approach. 

We also derive an universal mass formula for all Hirota's solutions to the Affine Toda 
model valid for all underlying Lie groups. Embedding of the Affine Toda model in the 
Conformal Affine Toda model plays a crucial role in this analysis. 
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1 Introduction 



Toda field theories are among the most important examples to study integrability and con- 
formal invariance. They can be classified according to an underlying Lie algebra. The 
conformal Toda (CT) theories are conformally invariant and completely integrable models 
associated to finite dimensional Lie algebra Q. The affine Toda (AT) models can be regarded 
as the conformal Toda perturbed in such a way that although conformal invariance is broken, 
complete integrability is preserved. The underlying Lie algebra is the Loop algebra Q and 
these models are known to possess soliton solutions. The third class, the conformal affine 
Toda (CAT) theories, are obtained by recovering the conformal invariance in the AT models 
by introducing two extra fields. These fields can be shown to be related to adding central 
extensions to the Loop algebra leading to a full Kac-Moody algebra underlying the CAT 
model. 

The interest in the soliton solutions of the Affine Toda (AT) field theoretical model is 
partly motivated by the well-known results for the simpler case of the Sine Gordon field 
theory and partly by the special status of the AT theory in view of its interpretation as 
deformed unitary minimal theory. 

This paper is devoted to obtaining, in a systematic manner, soliton solutions of the AT 
model and linking them to solutions and properties of the CAT. This is accomplished by 
using the Hirota's tau-function method which was already successfully been apphed to the 
AT field theory, although not all solutions were found, to the cases of the Lie algebras Ar 

Cr ^ and other algebras The method involves introducing one extra tau-function 
apart from those assigned to each AT field. It was within the context of the CAT model 
in that this fact was made clear in view of the new fields introduced in order to restore 
conformal invariance. 

There are three main (but related) goals which we have accomplished in this paper: (i) A 
description of soliton solutions obtained by Hirota's method for the AT model defined on an 
arbitrary Lie algebra, (ii) Analysis of the Hirota's method as a perturbation method based 
on recurrence relations, (in) Application of the link between Conformal Affine Toda (CAT) 
model and its integrable deformation represented by AT theory to reveal new features of the 
latter. One of the results we obtain here from general arguments of conformal field theory, is 
a general mass-formula valid for soliton solutions. A parallel and complementary approach 
to find soliton solutions for the AT model was developed in [Q using the Leznov-Saveliev 
method, see also |^ for the Baklund's method approach to the problem in the case of A^. 

The Hirota's method employed in this paper constitute essentially a perturbative ap- 
proach based on recurrence relations which can be conveniently characterized in terms of the 
eigenvalue problem of a matrix Lij = ifKij where is the extended Cartan matrix and 
if are integers in the expansion of the highest coroot in terms of simple coroots ^ of an 
algebra Q. A soliton solution is assigned to each eigenvalue of Lij and there can be at most 
rank Q + 1 when no degeneracy is present. The fact that the is singular is crucial in 
truncating the perturbative series leading to an exact solution. For the purpose of this study 
we have established a simple way of solving this eigenvalue problem for L^j matrix in terms of 
Chebyshev polynomials and their special properties. The situation is changed however when 
there is a degeneracy. There are two separate sources of degeneracy. One comes from the 
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degeneracy of the matrix Lij. Any linear combination of the degenerate eigenvectors gives 
rise to a sohton solution within the Hirota's method. This feature clearly demonstrates the 
nonlinear superposition principle underlying the method. Another degeneracy is an intrinsic 
feature of Hirota's perturbation scheme itself. The relevant recurrence relation is based on 
adding to the n-th order perturbed solution non trivial elements of the kernel of Lij — n^X6ij 
(for A an eigenvalue of L) which exists in some cases {SU{6p) and SP{3p) with p integer). 
Remarkably the Hirota series truncates in these cases too, producing a new class of solutions. 
It is generally true that the sohton solutions obtained in the degenerate cases truncate at 
higher orders of the tau-functions. 

There are two related ways of understanding AT model and its relation to CAT model. 
One way of thinking introduces AT model as a deformation of the CAT model preserv- 
ing integrability of the theory. Such a deformation modifies the underlying W^o symmetry 
structure of CAT but still maintain enough symmetry to keep the AT massive structure 
both attractive and tractable for physicists. On the other hand there is also a parallel ap- 
proach, introduced first into the literature in 0], which views AT model as a version of the 
CAT model with conformal symmetry being gauge-fixed. This second approach proves to 
be very useful in this paper defining the precise limit in which one model goes into another. 
One particular consequence of the connection between the CAT and AT models is a special 
simple relation between their respective Energy-Momentum (EM) tensors. The modified 
EM tensor of CAT model in some special limit consists of a sum of EM tensor of AT 
model and the extra surface term. Since there is no mass scale in the CAT model the only 
contribution to the soliton masses of the AT model comes from the pure surface term which 
ensures enormous simplicity of the final universal expression for the soliton masses. 

The outline of the paper is as follows. In Section 2 we discuss a connection between AT 
and CAT models and relation between their EM tensors. In Section 3 the Hirota's method 
is discussed as a perturbation approach and the general features of this perturbation and 
recurrence relations are described. We also derive the general mass formula and discuss 
topological charges. In Sections 4, 5, 6 and 7 this theory is used to find and describe soliton 
solutions of the AT model with the underlying Lie algebras Ar = SU{r + 1), Cr = Sp{r), 
Dr = S0{2r) and B,. = S0{2r + 1), respectively. In Sections 8, 9 and 10 this analysis 
is repeated for AT models with exceptional Lie algebras G21 -F4 and finally we make a few 
comments for the En with n = 6, 7, 8. In Appendix A we list for reader's convenience a 
number of technical properties of Chebyshev's polynomials. 



The equations of motion of the CAT model associated with a simple Lie algebra Q are given 



2 The CAT and AT Models 



by 0, S,0: 




(2.1) 
(2.2) 







(2.3) 
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where Kah = '^cia-Oihl is the Cartan Matrix of Q, a,b = 1, rank^ = r, ip is the highest 
root of G, = 2il).ai,/al, If are positive integers appearing in the expansion ^ = /^^, 
where are the simple roots of Q and g", and q are couphng constants. The derivatives 
d± are w.r.t. the hght cone coordinates = x ±t. 

These equations are invariant under the conformal transformations 



x+ ^ x+ = f{x+) , x_ x_ = g{x^) (2.4) 



and 



^-u(x+,x.) _^ ^-u(x+,x^) ^ ^_^-^B ^_^^B ^-v{x+,x.) ^2 6) 

....... . 

dxj^ dx 

where / and g are analytic functions and B is an arbitrary number. Therefore e'^" are scalars 
under conformal transformations and can also be arranged to be a scalar by setting 5 = 



The equations ( p.l|) -( p!3D can be written in the form of a zero curvature condition (for 
the associated linear system) [Q, ^ 



d+A^ - (9_v4+ + [v4+, AJ\ = (2.8) 



where 



and 



A+ = a+$ + e'"^^8+ , A_ = -a_<l> + e-'^'^*£_ (2.9) 



- /rankQ \ 



(2.10) 



rankQ rankQ 

^+ = E > E <l''E\+q'E--' (2.11) 

a=l a=l 

where T3 = 26. H'^ + /iD, with = | Z]a>o ^ (with a being the positive roots of G), is the 
generator used to perform the so called homogeneous grading of a Kac-Moody algebra. H^, 
D and C are the generators of the Cartan subalgebra of the affine Kac-Moody algebra Q 
associated to Q. E^^ {E^^J, a = 1, 2, . . ., rank^ = r, E^^ (E^-^) are the positive (negative) 

simple root step operators of Q. 

The Lagrangian for the CAT model is given by 

^ = T E ^Ka^d.^'^d"^' + -^d.^'-d^r^ + <f\d,^d^u - U{^, (2.12) 

^ a,b=l ^ a=\ ^a ^ 
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where h is the Coxeter number of Q and the potential being 



U{^, r/) = V l^e'^(^-^''+'^) + ^e^i-K^''^'+ri) (2.13) 
^1 ^ 



Let us introduce, for convenience the vector 



^ - E ^V'" (2.14) 



a=l «a 



The potential (|2.13|) can then be written as 



U{^, r]) = j2 ^e5>^-^+'') (2.15) 
where we have denoted = —ip as the extra simple root of the affine Kac-Moody algebra 

g. 

The potential ( p. 151) is invariant under the transformation 

if ^ if + ; rj^Tj + ^n (2.16) 

q q 

where ii" is a coweight of i.e. /x*" = Sa=i"^a^ where and n are integers. Such 
transformation is complex and therefore in the regime where the coupling constant q is pure 
imaginary the real vacuum solutions are degenerate. This generalizes to any simple Lie 
algebra the degenerate vacua of the sine Gordon model where the minima of the potential 
are identified with (co-) weight lattice of SU{2). Such degenerate vacua are responsible for 
the appearance of topological soliton solutions in the AT and CAT models as we describe 
below. 

Notice that the potential ( |2.15| ) and the CAT model equations of motion ( |2.1| )-( pl3| ) are 
invariant under the transformation 

+ a ; Ti^r]-C ; v^v + hC (2.17) 
where h is an arbitrary constant and C, is an harmonic function, i.e., = 0, and 



r 9), 

i^H-J (2-18) 



i=o 



where \j are the fundamental weights of the Kac-Moody algebra Q, i.e. Aq = (0, ^/q'O)' 
Aa = (Aa, \lt->^)-> a = 1,2, ...,r, with Aa being the fundamental weights of the simple Lie 
algebra Q and the scalar product defined as a.h = a.b + acbn + aobc for a = {a,ac,aD) 
and b = {b, be, bo)- When evaluating the scalar product we consider the simple roots aj as 
r + 2-component vectors, i.e. = {c(a, 0, 0), a = 1, 2, r and ao = {—ip, 0, 1). 



4 



The canonical energy-momentum tensor corresponding to the CAT model Lagrangian 
( p.l2| ) is given by 



^2 r 



a,b=l <^ 



1 



+ 9pa 



(2.19) 



where p, a = 0, 1 are space-time indices {do = dt, di = dx) and g^o = 1, gu = — 1, 
5^01 = 910 = 0. 

One finds that this E-M tensor has a non-zero trace given by 

e^ = 2f/(^,r^) (2.20) 

Since the CAT model is conformally invariant its E-M tensor can be made traceless. The 
usual procedure is to add a term Wpa^ = {dpd^ — gp^d"^) f with an arbitrary function /. 
Clearly the addition of such a term does not violate the conservation of Qpa- The trace 
of the extra term is given by Wj^ = —d'^j and this fixes the choice of / which renders the 
modified E-M tensor traceless. Accordingly one verifies that the modified traceless E-M 
tensor of the CAT model is given by 

= - ^ {^^^ - 9p<rd') f E + hu] (2.21) 

\a=l ) 

The CAT model can be obtained via a Hamiltonian reduction procedure from the two-loop 
WZNW model 0. The modification of the E-M tensor described above is equivalent to the 
one performed on the Sugawara tensor of the two-loop WZNW model in order for reduction 
procedure to respect the conformal invariance [Q. 

The Affine Toda model (AT) associated to a simple Lie algebra Q possesses rank-^ 
fields <y9", a = 1,2, ... rank-^; and it is not conformally invariant. Its equations of motions 
correspond to eq. ( p.l|) for r] = 0. As shown in 0] the AT models constitute a "gauge fixed" 
version of the CAT models. The t] field can be "gauged away" by a conformal transformation 
( p.4|) with /' = e''+*-^+-' and g' = e^~^^~^ where and rj- are the parameters of the solution 
of the free field rj, i.e. ti{x+, x_) = ■?7+(x+) + r]^{x-). Therefore for every regular solution of 
the r] field the CAT model defined on a space-time {x+,X-) corresponds to the AT model 
with the extra field u defined on a space-time (x+,x_) where 5+ = f'^'^ dy^e'^+^y+^ and 
x_ = dy-e^-^y-\ Such connection allows one to calculate exactly the perturbation of 
the 7] field on the dynamics of the AT model 0. 

For the particular solution 77 = the CAT and AT are defined on the same space-time 
and the dynamics of the fields if"" are the same on both models. In such case the E-M tensor 
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of these theories are related by 



= e^/^ |,=o +g [OA - 9p^d') ij2-^^ + hu] (2.22) 

\a=l / 

This relation provides a very interesting insight on the role of the broken conformal symmetry 
of the AT models. If one considers classical solutions of soliton type which can be put at rest 
at some Lorentz frame, then the energy of the solution can be interpreted as the rest mass of 
the soliton. Such mass should be proportional to some mass scale of the theory. Due to the 
conformal (scale) invariance such mass scale does not exist in the CAT model and therefore 
the soliton mass should vanish. Alternatively, this can be seen by taking the matrix element 
of Qp^^ between one soliton states of momentum p and p', obtaining p| 

{p\ Jdx e^/^ I p') = APpP^ + B [KpK, - K^gp„) (2.23) 



where P = and K = p — p' . Relation (|2.23| ) follows from the symmetry of Q^^'^ and 



its conservation d'^Q'^^^ = 0. If in the classical limit, which is obtained by considering the 
tree diagrams of the vertex ( 2.23| ), B has no poles at = 0, we shall have in the forward 



direction 



{p\ I dxQ';^ \ p) = Appp, (2.24) 



and because 6?^^ is traceless it follows we must have either = or A = for p^ ^ 0. 
Therefore does not give any contribution for the case of a soliton of mass different 

from zero {p^ ^ 0). 

However in the AT model the same is not true, and from ( p.22| ) one observes that the 



contribution to the soliton mass in the AT models comes from a total divergence contained 
in the second term of ( p^.22[ ). Indeed, denoting by M the soliton mass and v its velocity (in 



units of light velocity) one gets 

Mv /■°° 



AT 

ax c 

— oo 



dx 9oi 



roc. ( 2 \ 

q / dxd^dt ^'^^ + (2-25) 

\a=l ) 




\a=l "a 

This universal formula will be used to determine the masses of the solitons we will find below, 
using Hirota's method. This result was reported independently by Olive and we learned 
while typing this paper |^ that his results was obtained using very similar arguments to 
ours. 
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3 Hirota's Method 



We now construct solitons solutions for CAT and AT models associated to any simple Lie 
algebra Q using the Hirota's method |]10|- We introduce the r-functions as 

u=-^ia-\nTo) (3.1) 




E In ( ^ 1 (3.2) 



rankg /^Oli^' 
6=1 v7 / 

where R is the matrix with entries Rab = Sab+nf, with nf being the integers in the expansion 

ip = nfaa, and so = ^nf, and 

iR-\b = Sab-^ (3.3) 



Substituting these definitions in ( p.l| ) and (|2.3|) for r] = one obtains that the resulting 
equations can be decoupled into 

(rankQ \ 
1 - n r;''^'' (3-4) 
fc=0 / 

d+d-a = (3 (3.5) 

where Kij, i,j = 0, 1, 2, rank-^, is the extended Cartan Matrix of the Affine Kac-Moody 
algebra Q associated to Q (obtained from the ordinary Cartan matrix of Q by adding one 
extra row and column corresponding to the simple root ao = — ^). Furthermore 

Am.9,9_lnf^^^-M|z^ ,3,) 



and 

= 

I 



P = SLe^,A foranyj = 0,l,...,r (3.7) 



From (^j) and ( p.6[ ) one sees that, for B = 0, the r-functions are primary fields of conformal weight 

(0,0). Notice this definition of the r-functions is related to that of ref. by tj {^jf^ ■ Also, these 

r-functions are related to the Leznov-Saveliev jl^ solution of the CAT model Comparing (3J^) with eqs. 
(14) and (15) of ref. |^ (where g = 1) one can write (on shell and for ?/ = 0) 



(A(,) I e^+(-+)Af+(a;+)Afri(.T_)e-^-(--) | A(,)) 



To 



= e^"(A(o) I e^+(-+)Af+(a;+)Afr^(a;_)e-^-(---) | A(o)) 
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is a constant independent of the index j [Q. We have also set -do = 0. In the calculation we 
have used the fact that if , j = 0,1,2, rank-^, with /q = 1, constitute a null vector of the 
extended Cartan matrix. Indeed, ELo^ii^f = + Ea=i ^^If = 0, since ao = -^J 

and^ = EU/^ff 

The solution to ( p.5|) gives a as 

(r{x+,x^) = (3x+x^ + F{x+) +G{x^) (3.8) 

with F and G being arbitrary functions. 

Reference [|l| provided for the first time the Hirota's solution for the Affine Toda models 
associated to = SU{r + 1). The technical manipulations, required in order to write down 
the consistent Hirota's equation for the AT model, relied on the r-functions which exceeded 
by one the number of fields physically associated to the model. The origin of the extra r- 
function, namely tq, can be traced back, in view of the analysis in 0, to the u field revealing 
intrinsic connection of the Hirota's equations ( pj.4| ) to the structure of the CAT model. 

In the spirit of the Hirota's method we expand the r-functions in a formal power series 
in a parameter e as 

r. = l + 6..(^) + ... + e^-rP (3.9) 

We will be interested in solutions were the space-time dependence of the r-functions is 
given by 

with 

r = 7+X+ + 7_x_ + ^ = 7(x - t;t) + ^ (3.11) 

where (5-"^ are constant vectors to be found from Hirota's equations, and 7+ = ^{1 + v), 
7_ = |(1 — f ) and ^ are parameters of the solution. 

The basic idea is to expand Hirota's eqs. ( |3.4| ) in powers of e and solve them order by 
order. The method gives an exact solution if the series truncates at some finite order in e. 
The actual value of the parameter e is irrelevant in the procedure and it can be set to unity 
at the end of the calculation. Below we give a detailed procedure to solve Hirota's equations. 
But before that we would like to discuss some general properties of the method. 

Theorem 3.1 Let us suppose one has a solution of the Hirota's eqs. and let us denote by 
Ni the highest power of e in the Hirota's expansion of Ti, i = 0,1,2, ... , rank-Q. Then Ni 
are the components of a null vector of the extended Cartan matrix: 

KijNj = (3.12) 

Since the zero eigenvalue of is non degenerate and since if is a null vector of K, it follows 
that Ni = ntf , where k is some positive integer. 

Proof. Since Ka = 2 and Kij < for i ^ j one observes that by multiplying both sides 
of ( p. 41 ) by Tj the powers of r's will all be positive. After multiplication the term of highest 

power in e on the l.h.s. of (|3.4| ) will be e^*^" (rj^'^'^^ A (rj^'^''^^ which is zero because of the 
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ansatz ( ^.10 ). Notice now that on the r.h.s there is only one unique term contributing to 
the highest power of e leading to 



rankQ 

^ -Q (3.13) 

k=0,k^i 

Since the sum of the powers of e on the r.h.s. of the above equation should equal NiKa one 
obtains ( CT ) 

Therefore a necessary condition for the expansion ( |3.9| ) to truncate at some finite order 
(when ansatz ( p.lO|) is used) is that the matrix Kij appearing in the Hirota's eqs. (|3.4| ) 



should be singular. 

Substituting the expansion (|3.9| ) into Hirota's eqs. ( p.4|) and using ( p.lO|) and (|3.11|) one 



obtains that the resulting equation in order zero in e is trivially satisfied whilst the first order 
equation leads to 

Li/p = a4'^ (3.14) 

where Lij = if Kij, and A = -^^^^j^ = ^ ^]^^ ^ ■ Therefore the parameters of the solution are 
restricted by the possible eigenvalues of the matrix L^j. As we shall show below, except for 
the algebras SU{Qp) and SP{3p) with p a positive integer, the higher terms (J*-"^ {n > 2), are 
uniquely determined by 5^^\ Therefore, if the eigenvalues are non degenerate there can be at 
most rank-^ + 1 one-soliton solutions. However if there are degeneracies we can have many 
more solutions. Since a right null vector of Kij is also a right null vector of Lij and vice 
versa, it follows that there will always be a zero eigenvalue. This eigenvalue is not degenerate 
and the corresponding soliton solution is trivial in the sense that the cp^s fields (but not u) 
are constants. 



3.1 Soliton Masses 

From ( |3.1| ) one gets 

^ (t + hA=-t^ ^ogr, + ^ha + ± (3.15) 

\a=l "a / j=0 r a=l "a 

where we have used the fact that Y7a=i ~ ^ -'-)• Therefore from the relation ( p. 25 ) 
one obtains 



where we have chosen the functions F and G in the definition ( |3.8| ) in such way that the 
contribution to the mass from the a field vanishes. Substituting ( |3.9|) and ( 3.10 ) into ( p.l6 ) 
one gets 

v/T^ ^""h^^ l + 65f)er + ... + e^^5p)e^«r 

Recalling that F = 7(x — vt) + ^ we see that by taking 7 > the lower limit x — > — 00 does 
not contribute and the limit x — > cxd is finite. For 7 < the converse happens. Therefore the 
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mass is proportional to | 7 |. From Theorem ( |3.1| ) we have A'^- = kI^ where k is a positive 
integer. Therefore 

^HoM«E^iV,H.>^ (3.18) 

with h being the Coxeter number of Q. Recalling that A = 7^(1 — f^)/4/5 we therefore arrive 
at the following expression for the soliton masses in the AT model 

M = -^mVX 3.19 

where we have denoted m = y/j3. This mass formula is true for any algebra and for any 
solution constructed using the ansatz ( p.lOj ). Notice that the masses of solitons associated 
with a given eigenvalue A are quantized in units of ^my/^. In addition the soliton masses 
are proportional to the masses of the fundamental particles of the AT model. 

3.2 Soliton charges 

An important consequence of Theorem ( p.l|) is that the asymptotic values of the (f^s fields 
are always finite for the solutions constructed through the ansatz (|3.1CI|) . From ( p.l|) , (|3.9| ) 



and ( |3.1CI|) we see that, for 7 > 0, ipa{— 00) = y and since Ni = nlf with Zg = 1 we see that 



the limit of Ta/ro° as x — > 00 is finite and so from ([Ol) (pa{oo) is also finiteQ. For 7 < the 
limits are interchanged. Since the masses of the solitons, as shown above are finite, it follows 
that such asymptotic values of ip's are vacua of the potential ( |2.15|) . As discussed before, in 



the regime where q is purely imaginary the real minima of the potential are degenerate. We 
then introduce the topological charge of the solitons as {iq = q) 

Q = ^ r dxd^if = ^ (¥.(00) - </^(-oo)) (3.20) 

ZTT J-00 2,71 

The actual calculation of the charges presents some ambiguities which we now discuss. From 
(^, and we have 



2^ 1 20- r 1 20" 

where = ^ YZi=i ^^"^a- The r-functions are in general complex and therefore the logarithm 
function may not be single valued. One could therefore use the prescription, as in that 
log z = log I z I +i6 for z =\ z \ e*^, where < 9 < 2-71 and 9 — 9 = multiple of 2tt . However 
the implementation of such prescription requires also a criterion on the use of the identity 
log AB = log A + log B. The reason is that the result depends on the order in which the bar 
prescrition and the identity are used. For instance, one could get log AB = log A + log B = 
log I A I +log I B I +iaTgA + iaigB or log AB = log | AB \ +i(aigA + aigB). Due to 
such ambiguities the values of the charges one obtains will differ by a sum of co-roots of Q. 

We would like to point out that for x ±oo, Tq/tq" becomes a phase, i.e. lini2;_^±oo | Ta/TQ" \— 1 
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Therefore it is quite clear which coset in A|^/A^ the charge hes (where and A^j denote 
the CO- weight and co-root lattices of Q respectively), but the determination of the possible 
values of the charge, inside the coset, that a given sohton solution can have is quite delicate. 
For this reason we will give below the charges of the solutions we constructed, up to a sum 
of co-roots of Q. Following ( p.20| ) and (|3.21|) we will use the formula (for 7 > 0) 



Q 



1 2n T 

- I™ Z. — log — 



^0 



1 ^2a„, ^i^'') 
Z^ — log 



(3.22) 



where Nj is the highest power of e in the Hirota's expansion of tj, j = 0,1,2..., r. For 7 < 
the sign of the charge reverses. 



3.3 Recurrence Method 

We now describe a method to solve Hirota's equations recursively. We show how to determine 
the higher 5*^"^ (n > 2) from 5^^\ We write Hirota's equation ( p. 41) as: 



by introducing the auxiliary quantities: 



r A [tA 



^ I ^2 



and possessing the following e expansion: 



Gi 
F 



rankQ 

n 

k=0,k^i 



2r?(2) 



Clearly (|3.23|) must hold for each order i.e. 
where both sides of ( |3.26| ) can be rewritten as 



An) 



...,T 



(1) 



(3.23) 

(3.24) 

(3.25) 



(3.26) 

(3.27) 
(3.28) 



where Lij = ifKij and A^"''^ and S^"^ are sums of products of r's such that the sum of orders 
of r's for each term is n. From now on we choose the ansatz (|3.10|) for the r-functions. 
Substituting i ^27\ j and ([OSl ) into we get 



(n-l) 



n = 1,2. . . 



(3.29) 
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where A = 7-|.7_//5 and where we grouped together the lower order terms S^'^^^s with k < 
into one term V^"~^^ nonhnear in (5's and defined by 



n 



yin-1) ^ ^(l)^ ^ A, /fJ 15, ^ _ („„i) 

^ 'I exp 



For convenience we have introduced above new quantities = S'f'^ j expriF and ai = 

A-"V (7+7- exp nr). Notice that 1//""^=°) = due to: 

aP=bi'^ = ; z = 0,...,r (3.31) 

Hence the expansion of the r function must always start at the first order with 6^^^ being an 
eigenvector of Lij matrix as we explained in ( p. 141) . This forces the parameter A of Hirota's 
equation to be one of the eigenvalues A^^' of the Lij matrix. Since v}"^'^^ depends on d^''^ 
with A; < n as seen from its definition (|3.30|) , equation (|3.29|) is a recursive relation and can 
be used to determine higher order 5's. If n^X is not equal to any eigenvalue of Lij then the 
recurrence relation becomes 



4" ) = 5("\:;.V/'^-') ; Slf ^U,-n'X6ij (3.32) 

If n^A equals one of the eigenvalues of L^j we have to be more careful. In fact, such kind of 
"degeneracy" opens 
both and (5f 



"degeneracy" opens the way for constructing new solutions of Hirota's equation. Expanding 



4"^ = E</-f' ; vt-'^ = T.4^'^^?^ (3-33) 

k k 

in terms of the eigenvectors vf^^ of the L^j matrix: 

L.^vP = Xi%p (3.34) 
and plugging this expansion back into eq. ( p.29|) results in 

(aW - n^X) df:; = c[rf ^) (3.35) 

So, if A^^l 7^ n'^X then d\j^^ is determined from ( |3.35| ). However if A^^' = ra^A, is unde- 
termined and we can choose it arbitrarily. Notice however that for consistency Cj^]"^'' has to 

vanish in such case (in fact, in all cases where the kernel is non trivial we found that Cj^] 
consistenly vanishes). Such non trivial kernel of L^j — li^XSij introduces new parameters in 
the procedure and we will show below that it leads to extra soliton solutions of the AT and 
CAT models. 

As we are going to show in the following sections, the eigenvalues of the matrix Lij for 
the classical algebras can be written as 

A[^]=4cesin2(^^) ; j = 0,l,2,...,rg (3.36) 
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where cg = 1 and rg = rank-^ for Q = SU (r + 1), and Sp{r), and eg = 2 for ^ = S0{2r) and 
S0{2r + 1), whilst rg = rank-^ - 1 for ^ = S0{2r + 1) and rg = rank-^ - 2 for ^ = S0{2r). 
S0{2r + 1) has an eigenvalue A = 2 and S0{2r) has two eigenvalues A = 2. 

Consequently the kernel of Lij — n^Xdij for the classical algebras is non trivial when there 
exists j, f and n such that sin^ {^^^ = sin^ {,^\ '^^^ solutions for such equation are 
very scarce and can be obtained from the following more general theorem [T^ which we here 
present without proof. 

Theorem 3.2 The only solutions of the equation 

sin (avr) = csin (fevr) (3.37) 

with a, b and c non-zero rational numbers are when either c = ±1 or {sin (air) , sin (bn)} C 
II _i 1 _i| 

Therefore the kernel is non trivial for SU{Qp) and Sp{3p) only, with p a positive integer, 
(where in both cases h = Qp) with j = 3p, j' = p and n = 2. We will show later that in such 
cases we obtain new soliton solutions. 

For the exceptional algebras one can check that the kernel of Lij — n'^X6ij is always trivial. 

Summarizing, using the above method we can find the ^^"^ recursively using ( p.32| ) and/or 
( p.35| ). We expect to find the soliton solutions for the cases where this series truncates. That 



is, we have to have: 



y/""^^ = ioT n> Ni (3.38) 
since from now on all S^"'^ will be zero. 



3.4 Calculation of 

Kn-l) 



(n-1) 



The calculation of defined in ( p.30 ) can be easily performed using the Leibniz rule 



^^^^ — Y4=o t4 d""-^ ^^^^ ^^^^ arbitrary function H of the tau-functions 



we have if^"^ = d^H/de'^\e=o/n\. Let us first give a derivation of used to determine 

y/""^^ in (g). From (^ and (§1^ we have Gi = nd+d^n - d+Tid.n. Therefore 

n 

Gt^ = 7+7- E ((^ - If - Kn - /)) ^^-""e^ (3.39) 

Dividing by 7+7- exp (tiF) and subtracting the term ra^^j-"^ corresponding to dj^d-Ti we 
obtain 

n-l 



= ^ - 3n' ' oi2\,{i),{n-i) 
1=1 



' = ^ (n^ - 3nl + 2f) 5^5^^ (3.40) 
1=1 



The final form of depends through b\ on the Cartan matrix and varies therefore 

from algebra to algebra. However there are few generic terms which always appear in -F/"'^ 
As seen from ( p.25|) these terms are: (rf )'^"\ {Tk-iTk+if^\ (TiTj)^"^ and {TiTjTkY"'\ Therefore 
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we give here the results for generic quadratic and cubic terms in r-functions. Applying the 
same procedure as for the Gi function we arrive at 

n-l 

(r,r,)("Ve"^ = + + E "'^ (3-41) 

1=1 

and similarly 

n—l n—l 

1=1 1=1 

n-l n-l l-l 

+ E S^^Sf^ + E E st^Sf-^6t'^ (3.42) 

1=1 1=1 m=l 

Based on these results we will complete the calculation of V^^'^^ below for the relevant 
algebras. 



4 A, - SU{r + r 

The extended Cartan Matrix for SU{r + 1) is given as 



/ 



K 



and Hirota's equation 





V -1 

reads 

= /5 



-1 







-1 \ 


2 -1 










-1 2 


-1 










-1 


2 


-1 










-1 2 J 


- 'Tj+l^j- 


-0 


for 


J = 0,1,2, 



(4.1) 



(4.2) 



where, due to the periodicity of the extended Dynkin diagram, it is understood that Tj+^+i = 
Tj. The first order Hirota solution is obtained from the eigenvectors of Lij = tf Kij as in 
( P.14| ), namely 

Since if 



LijVj 



with 



1 for all i's, this yields the system of equations 

Vj-i — (2 — X)vj + Vj+i = for j = 0, 1, r 

Vj+r+l = Vj 



(4.3) 

(4.4) 
(4.5) 



Eq. ( [4.4| ) can be recognized as the recurrence relation for the Chebyshev polynomials with 2 — 
X = 2x (see appendix A). Therefore any linear combination of the Chebyshev polynomials of 
type I and II satisfy (|4.4| ). The relation ([4.5|) is equivalent, in fact, to the secular equation for 
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L and determines the eigenvalues. Using the trigonometric representation of the Chebyshev 
polynomials given in appendix A we get 



Xj = 4sm 

J Vr + 



y) forj = 0,l,2,...,r (4.6) 



By inspection Xj = Xr+i-j which clearly indicates that all eigenvalues, apart from Aq = 
and A(r+i)/2 = 4 (for even r + 1), are degenerate. The corresponding eigenvectors are 



/ 2iTilk\ r\,i / 2'n'ilk 



= ^^P^T^j ' ^2]fc = exp (^-^j / = l,2,...,L2j 
= (-1)'= for r + 1 even (4.7) 



where k = 0,1,2, ... ,r and 



(r — l)/2, if r is odd 
r/2, otherwise. 



Let us now apply our Hirota perturbation method to SU{r + 1). From (|3.3CI|) , (|3.4CI|) and 
( ^.41| ) we find a closed expression for V^/"-') to be 



1=1 



According to (|3.14|) , 5^^^ has to be an eigenvector of L. We start with the eigenvalues 
Ao = and Ar+i = 4 (for r + 1 even) which are not degenerate, and so we have df'' = 1 

and 6j^^ = (— !)■' respectively. Substituting these vectors into expression ([4.8|) for n = 2 
we see that vj^^ = 0. Consequently from ( p.32| ) we have 6^^ = 0, which leads to vj'^^ = 

(3) 

and 6j=0 and so on. Hence the Hirota perturbation series truncates at first order. For 
Aq = we then get tj = 1 + e^, which is a trivial solution in the sense that the yj's fields are 
constants (see (|3.1|) ). For Ar+i = 4 = '^''"J" (for r + 1 even) we get 



r, = l + (-l)V (4.9) 

According to ( p.l9| ) the mass of such solution is M = ^^^^^rri (r + 1 even). 

For the remaining set of (degenerate) eigenvalues we consider the following general linear 
combination: 



= yi exp j + y, exp (^-^ j / = 1, 2, . . . , [f] (4.10) 

where j = 0, 1,2, ...,r, i/i and y2 are arbitrary constants. The calculation based on 
gives: 



(2nUj\ ( 2ml J 
+ vo exp 



Vj'^ = -2ym ( 1 - cos ( ^ ) ) = -yiy2M^ - \i) / = 1, 2, . . . , [|] (4.11) 
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where we used expression ( [4.6| ) for the eigenvalues A/. Rewriting ( [4. Ill ) as V^^^ = ?/i|/2A/(4 — 
Ai)f and expanding 5*^^-* = d^^\^'^^ in terms of the eigenvectors from ( [4 .7] ) we get after 
substitution into (|3.35| ) 

6^'^=ym{l-\\)v^'^ (4.12) 

Substituting this into ( [4.8| ) we find V"*^^-* = and consequently S^^'^ = 0, which in turn leads 
to 1^(3) = and 5^^^ = 0. One can verify that l^^") = with n > A. Therefore 5^") = for 
n>3. 

To summarize, the general solution of Hirota's equations for Ar ~ SU{r + 1) is given by 
the tau-function 

r, = 1 + (^y, exp (^^^ + exp ) + yiy2 (l - ^A,) e'^ (4.13) 

with r given by ( p.ll| ), -^^j^ = A/ = 4sin^^^^ and / = 1, ...,[|]. The solutions given 
in ( [4.13| ) contain those obtained by Hollowood when either ?/i = or ?/2 = where the 
corresponding masses are given by eqns. ( p. 19] ) with k = 1, i.e. 

8(r + l)^^.^ Jt^ (4.14) 

However when considering both yi and ?/2 7^ , the solution truncates only in the second 
order in e yielding a solution with mass twice as large. We should point out that in such 
case we have a "two-soliton-like" solution surviving in the static limit. 

If we consider the topological charge given in ( |3.22| ) we find that when yi = the 
corresponding charge is 

^ 1 2aa f -2TTla\ ^ , , , r\ 

for some /?_ G A^. Taking now y2 = the charge is 

for some /?+ G A/j and / = 1, [|], cjj z = 1, ...r are the fundamental weights of SU{r + 1). 

The case of solution( |4.13|) when yi,y2 7^ 0, under such prescription correspond to charges 
lying in the root lattice A^. 



4.1 Second type of degeneracy for SU{6p) 

As noted in section |3.3| , for the eigenvalue Ap = 1 of SU{6p) the higher order ^'^"'^'s {n > 2) 
are not uniquely determined by S^^^ since the kernel of Lij — n^Xdij is not trivial for n = 2. 
We therefore can add to 5^"^^ a term proportional to the eigenvector corresponding to Asp = 4, 
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introducing a new parameter in the solution. The first order contribution to the tau-function 
corresponding to Ap = 1 is 



b) = yi exp 



/ 2-nipj \ 
V 6p ) 



+ 1/2 exp 



6p I ^i^-P(— j+^2-P(- — 



(4.17) 



where yi and 2/2 are arbitrary constants reflecting the ordinary degeneracy between eigen- 
values Xp and Xep-p. From ( [4.11| ) we get 



Vj''^ = -3ym 
and it is easy to see that the linear combination 



5(2) 



-V 



[A=0] 



+ zv 



[A; 



3pJ 



3l/ll/2 



(4.18) 



(4.19) 



for arbitrary z satisfles the second order recurrence relation: {L -2^\pl)6^^^ = V^^^ recalling 
that 2^Ap = Asp. Plugging now (|4.19| ) into (^78]) we obtain 



(2) 



(-1)^-+^. (5« + 5ji + 5f\) = -2z {y^u\^jf + +y2V^f) (4.20) 



(2) 

hence Vj is a linear combination of two eigenvectors corresponding to the degenerate 
eigenvalue A2p = 3. From recurrence relation (L — 3'^XpI)6^^^ = V^^'' we now flnd 6f = 



(4) 



z{—iy6j /3. One more step of recurrence procedure yields V^^^ = —zyiy2{—iy and 5j 
{—ly zyiy2/l2. At this point the recurrence procedure terminates and the flnal solution is 
given by 



( (^'^3\ ( 

= 1 + I yi exp ( — 1 + ?/2 exp ( — — 

Z ■ ( /^TTjA 

+ gi-l) ( 1/1 exp ( — 1 +1/2 exp 



r 1/1 1/2 



.2r 



e^^ + (-1) 



,zyxy2 4r 

e 

12 



(4.21) 



which of course reproduces ( [4.13| ) for 2; — » 0. The mass in such case is given by the mass 
formula (|3.19|) with k = 4, namely M = ^^5^. The topological charge in this case can be 
verifled to lie in the coset uj-^p + of A^y/A^j. 



5 Sp(r) Case 

The Cartan matrix is given by 



K 



( 2 

-1 






V 



-2 
2 
-1 








-1 
2 



.. 

.. 

-1 . . 

- 







2 

-2 



\ 



-1 

2 / 



(5.1) 
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/f = 1 for z = 0, 1, 2..., r, and yields the following Hirota's equations 

^oA(ro) = f3{r^-r^) 

tI a {to) = P{rl- Ta+iTa^^ for = 1, 2, . . . , r - 1 

r2A(r,) = /5(r2-r^i) (5.2) 

The corresponding eigenvalue equation LijVj = KijVj = \vi reads in components as 

{2-\)vq-2vi = (5.3) 
-Va-i + {2 - X)va - Va+i = for a = 1 , 2, . . . , r - 1 (5.4) 
-2w^_i + (2 - A)i;,. = (5.5) 

From (|]1|), has to be one of the eigenvalues oi L = K. Using the results of section |3.4| 
the function V^*-^^ from (p.30|) reads 

Vi'\5^'^) = i6iy-i5[y (5.6) 
Vl'\S('^) = (5«)^-52i5«i fora = l,2,...,r-l (5.7) 

v;(^)(5«) = (5«)^-(e\)^ (5.8) 

Our soliton construction is based on two results: 

1. The eigenvalues of the extended Cartan matrix are Xj = 4 sin^ (j7r/2r), j = 0, . . . ,r 

2. Vl^\6^^^) = (1 - a;2)(5j^))2, where x = {2 - X)/2 and = for n > 3. 

The proof of these results is based on the basic properties of Chebyshev polynomials 
||T3| , |T^ listed in Appendix A. Let us start on the proof by first introducing the variable 
x = (2 - A)/2 into (|^, (U) and and rewriting 

Vi=Ti{x)vo ;i = 0,...,r (5.9) 



Now equations ( ^.3|) and (|5.4| ) are trivially satisfied due to the basic recurrence relation (|A.l 
The remaining equation ( |5.5| ) becomes 

xTr{x) - Tr-l{x) = Q (5.10) 

This equation is not satisfied automatically by the Chebyshev's polynomials. It is in fact 
equivalent to imposing the secular equation for the SP{r) Cartan matrix. Using ( [A.2| ) we 
see that ( ^.10|) is equivalent to 

{I - x^)Tl{x) = Q (5.11) 

So the (r + 1) solutions to ( |5.1CI| ) are given by the extrema of Tr from ( |A.5| ). This corresponds 
to 

A.=4sin^(g) (5.12) 
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which proves the first observation. Next, substituting 



into 



we get 



Using ( |A.3| ) we get for (|5.7|) 

Vl'^ = {n - T.+ir„_i) = -(To - T,M'>f = (1 - x%5\, 

Using ( |5.1(J| ) we get 

V!^'^ = {T^ - Tl,) = (1 - x^)T^{5^^^f = (1 - x-mi),) 



(1)a2 



(5.13) 



(5.14) 



(5.15) 



where we used the fact that eigenvalues correspond to extremas for % (see ( |A.6|) ). Therefore 
( p.l3| ) and ( |5.14] ) hold for any value of x, while ( p.l5| ) is valid only when x corresponds to an 
eigenvalue of the Sp{r) Cartan Matrix. So the first part of point 2 is proved. We now prove 
the rest. The eigenvector corresponding to A = (x = 1) is v^^^ = 1, \/i as can be seen from 



( |A.?| ) or direct inspection. Therefore V^"^^ is proportional to v"^ °. From ( p.33|) and ( p.35|) 



we see that 5^^^ is also proportional to v 



A=0 



5(2) 



X] 



(1)^2 \A=0 



X ^ 1 



(especially all components are equal). Note, that 5^"^^ vanishes for A = (x 
indeterminacy in ( p.35| ) we don't get this result from ( ^.16|) 
From direct calculation we get 



T^(3) 

^0 


= (2- 






^(3) 


= (2- 




^(1) ^(2) x(2) 


v;(3) 


= (2- 


A)5(i)5(2) 





l,2,...,r-l 



(5.16) 
1), due to 



(5.17) 
(5.18) 
(5.19) 



After factorizing the equal components of ^'■^^ we are left with the eigenvalue equation for 
6^^^ and so V^^^ vanishes. From (|3.35| ) also 6^^^ must vanish. Hence V^'^^ will be the same as 
V^"^^ with 6^^^ replaced by Again since all components of 6^"^^ are equal V'^'^^ vanishes. 
Since all V*^"^ are quadratic in 5's and since the highest nonvanishing S is 5^"^^ it follows that 

From the above considerations we can write the general soliton solution for SP{r) as 



'^l±^{6iyT,{l)e'^ x ^ 1 z = 0, 1, . . . , r (5.20) 



r, = 1 + T,(x)#^e^ + 



where x = (2 — A)/2 and A = 7+7_//5. Notices that for x = —1 the series breaks at the first 
order and this corresponds to known solutions of Sp{r) For x = 1 we have 



r, = l + r,(l)5«er 
and so all the if fields are constant since all r's are equal. 



(5.21) 
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Masses can be easily analyzed from the universal soliton mass formula ( 3.19| ). Clearly, 
K = 2 for the new solutions in (|5.20|) and taking into account the value of the Coxeter number 
for Sp{r) {h = 2r) we find 



16r / irj 
— — m sm — 

7/'2 \2r 



J = 0, 



(5.22) 



For solutions with a; = — 1 the mass would be half as large. 

The topological charge defined in ( |3.22| ) lies in the coroot lattice for x ^ —1 since 
all r's possess in common the same coefficient of the highest power in e^. For x = —1 the 
corresponding topological charges given from ( ^.20|) is 



1 2aa 
27r ai 



(5.23) 



where uji is a fundamental weight of Sp{r) leading to the defining representation. Again the 
solutions ( |5.2CI| ) provide solitons with topological charges lying in the coset A'^/A]^ where 
A^ and A^ denote the coweight and coroot lattices of Sp{r) respectively. 



6 Dr = S0{2r) 

The case of D4 = SO (8) has a special interest and so we start by discussing it separately. 



6.1 L»4 = S0{8) 

The matrix Lij = if Kij is given by 



/ 2 -1 \ 

2-100 

-2 -2 4 -2 -2 

0-120 

V -1 2 y 



(6.1) 



and the integers if are (1, 1, 2, 1, 1). The eigenvalues of the matrix L are (0, 2, 2, 2, 6) and 
the corresponding right eigenvectors are 



,A=0 



„A=2 



/ 1 \ 
1 

2 
1 

V 1 / 

The Hirota's equations read 



/ 


1 \ 




/ 


1 \ 






1 \ 




( 


1 \ 




-1 






1 






-1 






1 







■ 7,^=2 - 







■ 71^=2 - 











-4 




1 




-1 




-1 






1 


V 


-1 ) 










I 


1 ) 




\ 





(6.2) 



rl A (ro) 
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rl A (n) 
r| A {r,) 
rl A (ra) 
A (r4) 

The solution corresponding to = 



T2 

- T0T1T3T4 
^2 

r2 



(6.3) 



A = 6 is given 

= (1,1,-4,1,1) 
= (0,0,1,0,0) 



(6.4) 



and all the remaining (5's vanish. Therefore according to Theorem p.l| we have k = 1. So, 
from ( |3.19|) we get the mass of such soliton is M = ^Vl\f^m. 

The eigenvalue A = 2 has multiplicity 3. We then apply Hirota's method by taking 5^^^ 
to be a general linear combination of the three eigevectors. The solution is then given by 



^(1) 



-X\ + X2 - X3 



X\ — X2 — X3 
y -Xi - X2 + X3 ) 



/ X\X2^ X\ X3 + X2X3 \ 
-Xi X2 + X1X3- X2 X3 

3{Xi^ +X2'^ + X3^) 
-Xi X2 - XiXs + X2 X3 
Xi X2 — Xi X3 — X2 x^ 



= A 



/ 1 \ 

1 

-16 
1 
1 



( \ 


1 





/ \ 






p2 



( \ 


1 





(6.5) 



xx X2 + xx xz + X2 a;3^) /9. 



and the higher (5's vanish and where P3 = x\ X2 x^/27 and P4 
The solution holds true for any value (even complex) of the parameters Xi, X2 and X3 used 
in the linear combination. 

The properties of the solution seem to be insensitive to the actual values of the parameters 
x's as long as they do not vanish. However when one or two of them vanish the solution 
changes substantially. According to Theorem p.l| we have k, = 1,2,3 when two, one and 
none of the x's vanish respectively. Then from (|3.19|) the masses are given by 



Ml 
M2 
M3 



^12y2m 
^24v^ 
^36v^ 



m 



m 



when two x's vanish 

when one of the x's vanishes 

when none of the x's vanishes 



The topological charge defined in ( ^.22| ) yields for the solution ( |6.4| ), Q G A/j, the root 
lattice of SO {8). For those solutions given in (|6.5|), the topological charge depends upon the 
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parameters Xi,X2 and X3. By explicit calculation it can be shown that 

g(a;i,0,0) G Xs + Ar 

Q(0,X2,0) G X^ + Ar 

Q(0,0,X3) G Xs + Ar 

Q{xi,X2,X3) G Aij (6.6) 

where A^,, and A^ are the vector and the two spinor weights of 5*0(8) respectively. 



6.2 Dr = S0{2r) for r > 5 
The L matrix is here given by 



V 



2 

-2 










2 
-2 










-1 
-1 
4 
-2 





-2 

4 

4 
-2 








-2 

-2 
4 
-1 
-1 









-2 
2 




\ 






-2 

2 / 



(6.7) 



with the integers tf given by (1, 1, 2, 2, 2, 1, 1). To analyze the eigenvalue equation LijVj = 
Xvi we follow the usual procedure of first introducing variable x = (4 — A)/4 and expressing 
the components Vi in terms of Chebyshev polynomials. As usual the generic part of the 
eigenvalue problem takes the form of the recurrence relations for Chebyshev polynomials 
( |A.1| ). Trying to fit an ansatz in form of linear combination of Chebyshev polynomials into 
all eigenvalue equations produces the following result (in case A 7^ 2 or a; 7^ |): 

^'1 = ^^0 Va = 2 {Ua-1 - Ua-2) vq ] = 2, 3, . . . , r - 2 

Vr = Vr-1 = {Ur-2 " t^r-s) ^^0 (6.8) 

with the consistency condition, playing the role of the secular equation, 

Vr^2 = (4a; - 2) Vr or {x - l)T'._i{x) = (6.9) 
to be imposed on the solution (|6.8|). The solutions to ( |6l9| ) take the following form 



X, 



) sm 



A; = 0,l,...,r-2 (A ^ 2) 



2{r - 1)^ 

The case A = 2 has to be treated separately. One finds that the eigenvalue A 
degenerated with eigenvectors: 

/ 1 \ / \ 

-1 




(6.10) 
2 is twofold 



„,A=2 
^[1] 



A=2 

\2] 



\ 0/ 





-1 
V 1/ 



(6.11) 
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In case where r — 1 is a multiple of 3 (r — 1 = 3A^) degeneracy becomes threefold and the 
extra eigenvector has components given in terms of three components Vq, Vi, Vr-i playing role 
of free parameters 

V3+3k = V4+3k = -V6+3k = -V7+3k = -{vq + Vi) k = 0, 2, 4, 6, . . . 
V2+3k = k = 0,l,,---,N 

Vr = -t;._i + (-l)^+^(i;o + t^i) (6.12) 
The Hirota's equations read 



r'o A (ro) 


= P{rl- 




r'l A in) 


= P{r'i- 




rl A {r,) 


= 2p{rl 




rl A (rj 


= 2/?(rf 




A (r,_2) 


= 2P{tI 


2 ~ Tr-S^r-lTr) 


T^, A (r,,_i) 


= f3{rli 


- Tr-2j 


r'r A (r.) 


= P{rr 


- T-r-2) 



(6.13) 

Using the results of section |3.4| one obtains that the quantities y/"'^^'* defined in ( p.30| ) are 
given by 

.(n-l) _ v^/, , / 2 Q^; , o72n\ 



V^-'^ = -^(l-A(n2-3n/ + 2/2))5«5(" 



«=i 

n-1 

^^(n-1) ^ _^(l_^(^2_3^^^2/2))5P4"-') 
«=1 

n-1 

^("-1) = -5](2-A(n2-3n/ + 2/2))5«5?-') 
«=i 



n-1 / l-l ^ 

Z=l V m=l y 



T/(«-l) 
a 


n-1 

«=1 


- A(n^ 


' - 3n/ + 2/2 


)) ^i'^e-'^ 


T/("-l) 


n-1 

= -E(2- 

«=i 


-A(n2 


- 3nZ + 2/2) 


) a.e^'^ 




n-1 / 

+ 2EK2 

Z=l \ 


dn-l) 
3"r-l 






T/(«-l) 
*^r-l 


n-1 

«=i 


-A(r22 


- 3nZ + 2/2) 


) 5£ifc'^ 


T/{n-l) 


n-1 
Z=l 


-A(n2 


- 3n/ + 2/2) 





l-l 



m=l 



(6.14) 
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Applying the recurrence method described in section one can now construct the solutions. 

2{r-l) 



For the eigenvalues = 8 sin^ tT^ztt ~ "^^^s" 7^ 2, /c = 1, 2, r — 2, we obtain the solution 



To 



1 + 
1 , „r 



cos (2a — l)6k r or 
Ta = 1 + 2 ^ -^e^ + e^^ ; a = 2, 3, r - 2 

COS^fc 

Tr-l = l + {-l)'e^ 

Tr = l + (-l)'=e^ (6.15) 

where 9k = 2{r-i) where we have used the fact that the Chebyshev polynomials satisfy 

Ua-i{xk) — Ua-2{xk) = '^"'^iosg^^^'' with Xfc = cos2^^fc. Accordiug to Theorem |3]1| we see that 
K = 1 in this case and therefore form ( p.l9| ) the masses are given by 

2 kn 

Mk = — 8(r -l)V2m sin ; = 1, 2, r - 2 (6.16) 

ip"^ 2[r — 1) 

These solutions correspond to topological charges lying in the coset Qk ^ K + for k odd 
and Qk € for k even,where A^^ denote the vector weight of S0{2r). 

For A = 2 we have to consider the case when r — 1 is a multiple of 3 separately since 
as shown above the eigenvalue has multiplicity three. Let us consider first the case r — 1 ^ 
multiple of 3. The eigenvalue A = 2 has multiplicity 2 and the eigenvectors are given by 
( |6.11| ). So, taking 6^^^ to be a linear combination, with parameters yi and 1/2, of those 
eigenvectors and applying the procedure of section ^]3| we obtain the following solution 

To = 1 + yie^ + c{yi,y2)e'^^ 

n = 1 -yie^ + c{yuy2)e'^^ 

Ta = l + da{yi,y2)e^^ + c{yi,y2Ye'^^ ; a = 2,3,...,r-2 

r,_i = l-y2e^ + i-lY-'ciy,,y2)e^^ 

Tr = l+y2e^ + (-l)^-^c(yi,y2)e'^ (6.17) 

where F is defined in (|3ll ), A = 2 = c{yi,y2) = {yf + (-l)^'-^yi)/4(r - 1) and 

da = i-ir ((4(r - a) - 2) yf + (-l)'-(4a - 2)yl) /4(r - 1), a = 2, 3, r - 2. 

There are two special solutions which correspond to the choices y2 = ^i'"yi since in such 
case c{yi, 1/2) = 0. Then according to Theorem ^]T] we have k = 1 and from (|3.19|) the masses 
of these two solutions are the same and equal to 

2 r- 

Mi = — 4(r - l)V2m ; for y2 = ±i'yi (6.18) 
For y2 7^ ^i^yi we have k = 2 and therefore the mass is 



M2 = ^8(r - l)V2m ; for y2 ^ ±t'yi (6.19) 
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For y2 = fyi we can show that tha topological charge for either r = 2/c or r = 2A; + 1 lie in 
the cosets, Q € + A/j for even k oi Q E Xg + for k odd. When y2 = —i^yi the situation 
is reversed, i.e. Q e As + Ar for odd A; or Q G A^ + for k even. 

For the case where r — 1 = multiple of 3, the eigenvalue A = 2 has multiplicity 3. 
Analogously, we then take 6^^^ to be a linear combination of the three eigenvectors ( |6.11j ) and 
( |6.12| ) and apply the method described in section |3.3| . The calculations are a bit cumbersome 
and we give here the results for the simplest example corresponding to 5*0(14). The main 
properties of the solution for any S0{6p + 2) (pa positive integer) are the same for any 
p. Obviously as described above, 5*0(8) has also the eigenvalue A = 2 with multiplicity 
3. However for 50(8) the degeneracy of A = 2 is related to the symmetries of the Dynkin 
diagram whilst for 50(6j9 + 2) {p > 1) the degeneracy is accidental. Indeed, the solutions 
are different. 



6.3 50(14) solutions for A = 2 

The three eigenvectors for the eigenvalue A = 2 are 

(1,-1,0,0,0,0,0,0) 



,A=2 



.A=2 



(0,0,0,0,0,0,-1,1) 

(1,1,0,-2,-2,0,1,1) (6.20) 
By taking 5^^^ as a linear combination of these eigenvectors and applying the Hirota's pro- 



„A=2 



cedure described in section we get the following solution 



5(1) 



/ yi + y3 \ 
-yi + ys 



-2y3 
-2y3 


-y2 + z/3 

V 1/2 + Z/3 / 



/ ml -\_ ml -\_ yi y-i \ 

24 24 3 * 

ml _i_ ml. — ys 

24 ~'~ 24 3 
3^2 ,„2 



5(2) 



-7yi^ 

5yil 

-yi 



■\ + y3' 

5y2l 

13, 
7y2^ 



+ 



ys^ 

12 +y^' 

^ + 1/3^ 



4, ' ,4 

ml _i_ ml _ ?/2 ys 

24, ^ 24, 3 



5(4) 



\ 24 ^ 24 ^ 



5(3) 



^3 

27 



yi^+y2^ 



-2(1/1^+1/2^) 
23t;i^-13i/2^ 
92 

-13yi^+23t;2^ 

-2{y^'+y2') 



\ 



yi^+y2^ 



V 

















p, 1 yi'^ya'^ 
TD Tyi y-i 

^ , W 2 


+ 


. ^ 2 
5^2^ 2/3 

. W 2 


^4 + 




7 J/2 J/3 


D 2/1 

-^4 36 


+ 


?/2 3/3 
12 



















( ° ^ 




( ° \ 



















1 




1 
1 







1 
1 
1 



[0 J 



(6.21) 



where P4 = (yf + ylf /576, P5 = -?/3 (2/? + 2/1)' /2592 and = (^3 (2/? + vl) /216)'. 
We should distinguish four type of solutions: 
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1. We get three solutions by taking 1/3 = and 1/2 = ±iyi or 1/3 7^ and 1/1=2/2 = 0. 
According to Theorem ^]T] we have n = 1 and then from ( p.l9| ) the masses are 

Ml = -^24^/2771 (6.22) 

The topological charges for 7/2 = Wi can be shown to he in the coset Q E Xg + 
whilst for 1/2 = -iyi, Q e Xs + Ar. When ^ 0, and 2/1 = 2/2 = 0, Q G Ar. 

2. Taking 2/3 = and y2 7^ ±^2/i we have k, = 2 and then 

M2 = ^48^2 m (6.23) 

3. Taking 2/3 7^ and ^2 = ±^2/i again we have k, = 2 and then 

M3 = ^48V2m (6.24) 

4. Finally when 2/3 7^ and 1/2 7^ ±i|/iwe have k = 3 and then 

M4 = ^72V2m (6.25) 

In the cases 2, 3 and 4 above the topological charges lie in A^^. 



7 = SO(2r + 1 



7.1 50(7) 

The solitons solutions for the case = S0{7) are very similar to those of 5*0(8). This 
is a consequence of the fact that the Dynkin diagram of 5*0(7) can be obtained from that 
of 5*0(8) by a folding procedure. The Hirota's equations can be obtained from ( |0| ) by 
making T4 = and ignoring the equation for r4. The eigenvalues of the matrix L for 5*0(7) 
are (0,2,2,6). The solution for A = 6 for 5*0(7) is obtained from ( |6.4|) by deleting the last 
component of the vectors S^^^ and 6^"^^ Since the integers if are (1, 1, 2, 1) we get, according 
to theorem B.U that k = 1 for this solution. Then from ( p.l9| ) the mass of the soliton is 



The eigenvalue A = 2 of 5*0(7) has multiplicity 2 and the solution is constructed by taking 
5^^^ as a linear combination of the corresponding eigenvectors and applying the procedures 
of section p.3| . The solution one gets can be obtained from ( |6.5D by making 0:3 = xi and by 



ignoring the last component of all 5's (the solution will then depend upon two parameters 
xi and X2)- We have k = 1, 2, 3 for xi = and 0:2 7^ 0, xi 7^ and X2 = and a;i , 0:2 7^ 
respectively. The masses are then 



Ml 
M2 
M3 



^12V2m 
^24y2m 
^36\/2m 



when Xi = 0, 0:2 7^ 
when Xi 7^ 0, X2 = 
when xi , X2 7^ 



26 



7.2 Br 

The Lji = 



= S0{2r + 1) for r > 4 
ifKij matrix is here given by 



V 



2 

-2 









2 
-2 








-1 
-1 
4 
-2 






-2 
4 

-2 




-2 

4 
-2 








-2 
4 

0-1 



\ 







-4 

2/ 



(7.i: 



where the integers if are (1, 1, 2, 2, 2, 1). To analyze the eigenvalue equation LijVj = Xvi 
we follow the usual procedure of first introducing variable x = (4 — A)/4 and expressing 
the components Vi in terms of Chebyshev polynomials. As usual the generic part of the 
eigenvalue problem takes the form of the recurrence relations for Chebyshev polynomials 
( |A.1| ). Trying to fit an ansatz in form of linear combination of Chebyshev polynomials into 
all eignvalue equations produces the following result (in case X 2 or x 



Vi = Vo 

Va = 2 {Ua-1 - Ua-2) vq ; a = 2,3, . . . ,r - 1 

Vr = {Ur-1 - Ur-2) Vo 

with the consistency condition, playing the role of the secular equation, 

Vr-i = (4x — 2) Vr or (x — l)T/(x) = 
to be imposed on the solution (|7.2| ). The solutions to (|7.3| ) takes the following form 



A. 



I sm 




0,1, 



1 (A ^2) 



The case A = 2 has to be treated separately. One finds that the eigenvalue A 
corresponding eigenvector: 

/ 1 \ 

-1 



(7.2) 
(7.3) 

(7.4) 
2 has 



„A=2 

^[1] 







(7.5) 



V 0/ 

In case where r = 3(A^ + 1) with some = 1,2,... we get degeneracy with the extra 
eigenvector having components given in terms of two components vq, Vi playing a role of free 
parameters 



V3+3k 
V2+3k 

Vr 



Vi+Sk — 

k 



-V6+3k 

0,1,, 



-V7+3k 

.N 



■{vo + Vi] 



0,2,4,6, 



^ ' 2 



(7.6) 
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The Hirota's equations read 



2(3 [rj - roTiTa 



r„^A(rJ = 2/? r,^-r,_ir,+i a = 3, 4, r - 2 



rtiA(r,,_i) = 2/3 <_i-r,_2r, 



Tr-l 



(7.7) 



Using the results of section |3]^ one obtains that the quantities V^"' defined in ( p.30|) are 
given by 



Vr 



(n^l) 



n-1 



(n-1) 



(n-1) 



- ^ (l - X{n^ - 3nl + 2/2)j S^'S, 
1=1 

n-1 

-5:(l-A(n2-3n/ + 2/2))5«4"-') 
«=i 

n-1 

- ^ (2 - A(n2 - 3n/ + 2/^)) ^W^^""') 



«=i 

n-1 



i-1 



yin-i) 



T/(«-l) 
Vr-1 



/=1 
n-1 



m=l 



E ((2 - X{n^ - 3nl + 2f)) 5^^5^r'^ - 25ti& 
1=1 

n-1 

^ (2 - \{n^ - 3nl + 2f)) 5^1/^-1^ 



a = 3,4,...,r-2 



1=1 

n-1 



l-l 



1=1 

n-1 



m=l 



n-/) 



(7i 



«=i 



Applying the recurrence method described in section 3.3 one can now construct the solutions. 



For the eigenvalues = 8 sin 



2 k-K 7+7- 



2r /3 



7^ 2, = 1, 2, r — 1, we obtain the solution 



To = 1 + e^ 
Ti = l + e^ 

= 1 + 2^^^i^^^e^ + e^^ a = 2, 3, r - 1 

COSfc'fe 



(7.9) 



where (^k = ^ and where we have used the fact that the Chebyshev polynomials satisfy 
Ua-i{xk) — Ua-2{xk) = ^cos()l^^'° ^i^^ Xfc = cos2^fc. Accordiug to Theorem |3]l| we see that 
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K = 1 in this case and therefore form ( ^.19| ) the masses are given by 

— 8rv2m sin — 
ip'^ 2r 



Mfc = — 8rV2m sin— ; A; = 1, 2, r - 1 (7.10) 



and topological charges can be shown to lie in Q G A/j. Notice this solution can be obtained 
from ( |6.15| ) for S0{2r + 2) by a folding procedure. Identifying r^+i = Tr one observes that 
( |6.15| ) (for r —> r + 1) reduces to (|7.9| ) and the Hirota's equations (|6.13|) reduces to (|7.7|) . 

Consider now the case where r 7^ multiple of 3. The eigenvalue A = 2 is not degenerate 
in this case and taking S''^^ to be the eigenvector (|7.5|) and applying the procedure of section 
one obtains the following solution 

To = 1 + + 

4r 

r. = l + (-l)'^(l-^)e^^ + ^e^^ a = 2,3,...,r-l 
rr = 1 + ^-" (7-11) 



M = —8rV2m (7.12) 



where F is given by ( ^.11| ) with = 2. According to theorem |HT^ we have k = 2 and 
therefore from (|3.19|) the mass of such soliton is 

2 

and its topological charge lie in Q G Ar. Again, the above solution can be obtained from 
( |6.17| ) for S'0(2r + 2) (i.e. r — >• r + 1) by making yi = I and 2/2 = and identifying r^+i = r^. 

For the case where r = multiple of 3, the eigenvalue A = 2 has multiplicity 2. The 
solutions are obtained by taking 5^^^ to be a linear combination of the eigenvectors ( |7.5| ) 



and (fr6|) and applying the procedures of section |3]^. Like in the S0{6p + 2) case the 



calculations are cumbersome. One can in fact obtain such solutions for S0{6p + 1) from 
those of S0{6p + 2) by a folding procedure. For instance, the solution for 5*0(13) can 
be obtained from ( |6.21|) for the case of 5*0(14) by making 1/2 = and neglecting the last 
component of all S^^'^^s, n = 1, 2, 6. We then obtain three types of solutions: a) for yi = 
and ?/3 7^ we have k = 1 and the mass is Mi = :^24:^/2m■, b) for yi ^ and 1/3 = we 

have K = 2 and the mass is M2 = ^A'S\/2m and finally c) for j/i , j/3 7^ we have k = 3 and 

the mass is M3 = ^12\f2m. 

8 G2 

The matrix Lij = ifKy in this case is given by 

L = \ -2 4-6 1 (8.1^ 
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where the integers if are (1,2, 1). The eigenvalues of L are (0,2,6) and the corresponding 
eigenvectors are: 



,A=0 



The Hirota's equations for G2 are 

A (ro) 
A (n) 



2(3 (r^ - ToT^ 



For A = 6 the solution obtained, using ansatz ( p.lO|) and the procedure of section 
given by 



5(1) 




5(2) 




We therefore have k = \ and according to ( ^.19 ) the mass is 

2 



Ml 



ij)'- 



m 



For A = 2 the solution is given by 







/ 












V 










1/3 






^ 0, 






For such solution we have k = 3 and therefore from ( |3.19|) the mass is 



M2 = — 36y2m 



(8.2) 



(8.3) 



IS 



•4) 



.5) 



(8.6) 



Notice that these two solutions can be obtained from the 5*0(8) solutions by a folding 
procedure. Indeed, by identifying ri = rs = one observes that the Hirota's equations 
( |6.3| ) become (|8.3|) (after the relabeling ri T2). In addition, under such identifications the 

and the solution (|675| ) becomes (|8.6|) by seting xi = X2 = = 1. 



solution (|6.4|) becomes 
Both solutions yields topological charges lying in the co-root lattice of G2- 
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9 F4 

The L matrix for F4 is given by 



/ 2 
-2 



V 



-1 

4 -2 

-3 6 

-2 








-6 
4 
-1 



\ 



-2 

2/ 



(9i 



The integers If are (1,2,3,2,1). The eigenvalues of L are (0,6 + 2^3,6-273,3 + 73,3-73) 
and none of them are degenerate. 
The Hirota's equations read 



< A (ro) 


= f^{rl- 




rl A (rO 


= W[rl 




r^^ A (r^) 


= 3/3 (r^ 




r| A (rs) 


= 2/?(r| 


- ^2X4 


A (r4) 







(9.9) 

We construct the solutions using the ansatz ( |3.10|) and the procedure of section ^73| . The 
results are listed below. 

The solution for A = 6 + 2^3 = ^^^"^^ is given by 



^(1) 
5(2) 
5(3) 



(1, -2(2 + 73), 3(3 + 273), -2(2 + 73), 1 
(0,1,3(3 + 273) ,1,0 
(0,0,1,0,0) 



We have k = 1 and therefore from ( |3.19| ) the mass of the soliton is 

2 



Ml 



24V6 + 273 



m 



(9.10) 



(9.11) 



The solution for A = 6 — 273 



7+7- 



is 



5(1) : 

5(3) = (0,0,1,0,0) 
Again we have k = 1 and therefore the mass is 



(1, -2(2 - 73), 3(3 - 273), -2(2 - 73), 1 
0,1,3(3-273) ,1,0 



— r24V6 - 273m 



(9.12) 



(9.13) 
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The solution for A = 3 + -\/3 = jg 



where 



1^ 2 ' ' ' 2' 4 
5^') = (^ai,2 (27+ 14^3) ai,^,2 (3 + 2 Vs) ai,ai 

5^=^) = 6i(0,-2,8 (l + V3),l,0 

5^^) = a2(0,l,(63 + 4.3t),l,0 

^(5) = (0,0,0,0,0) 

5^^) = a? (0,0, 1,0,0) 



5 + 3i 

ai 



64 (71 + 41 V3 
(12417 + 7169 



128 (172947 + 99851 ^3) 
In this case we have k = 2 and therefore from (R.19I) the mass of the sohton 



M3 = -^48\/3 + VSm 



The solution for A = 3 — 



7+7- 

13 



IS 



2 ' ' ' 2' 4 



5(2) = f a2, 2 f27 - 14 VS) aa, |-, 2 f3 - 2 Vs) aa, 



64 

5^') = 62(0,-2,8 (1- Vs) ,1,0 

= a2(0,l,(63-4.3t),l,0 

5(5) = (0,0,0,0,0) 
5(6) = (0,0, 1,0,0) 



where 



-5 + 33 

02 = 



64 (-71 + 41 Vs) 

12417- 7169^3 
128 (-172947 + 99851 ^3' 
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Again k = 2 and the mass is 

M4 = -^AsJs- VSm (9.21) 
All solutions correspond to topological charges lying in the co-root lattice of F4. 

10 The En series 

The one-soliton solutions for the exceptional algebras Eq, E-j and E^ were presented in [|]. 
However not all static soliton solutions were obtained there. In the case of E^ the eigenvalues 

Ai = 6 — 2v^ and A2 = 6 + 2-\/3 have multiplicity 2. Therefore using our procedure of taking 

5(1) 

linear combination of the eigenvectors associated to the degenerate eigenvalue, one 
obtains new solutions for Eq. We do not present them here, but under a folding procedure 
they lead to the solitons of F4, constructed above, corresponding to the eigenvalues 6 — 2 a/3 
and 6 + 2^3. The masses of these new soliton solutions of Eq are the same as those of the 
corresponding solitons of F4. 



11 Conclusions 

The Hirota tau-function method employed in this paper has proved to be extremely powerful 
in providing explicit soliton solutions for the Toda models. In particular, when dealing with 
degenerate eigenvalues the method prescribes how to consider a superposition of solutions 
for the non linear problem. This yields, for instance, new solutions not previously discussed 
in the literature surviving in the static limit. 

We have also given the massess of such solutions by using a mass formula obtained by 
general arguments of conformal field theory. For many of our solutions we have discussed the 
topological charges and verified that they can be classified according to the cosets A^/A^, 
the quotient of the coweight and coroot lattices. 

It is still unclear, however, how to use the representation theory of Lie algebras in order 
to classify solutions and their topological charges as conjectured by Hollowood in [jT|. A 
very interesting and natural continuation of this work would be the study of the quantum 
theory of these solitons. This could be approached from several directions like bosonization, 
equivalence theorems and S-matrix p^ . 



A Chebyshev Polynomials. Definitions and Identities 

Here we list for completeness the basic properties of Chebyshev polynomials |]13|, |T^. They 
play an instrumental role in solving the eigenvalue problem of the Lij matrix for various Lie 
groups. 

The defining relation for Chebyshev polynomials of both types, T„(x) (type I) and 
Un{x) = T^^i/{n + 1) (type II) is given in terms of recurrence relations: 

Tn+i{x) - 2xTn{x) + Tn^i{x) = ] n>l 

Un+l{x)-2xUn{x)+Un^l{x) = (A.l) 
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with Tq{x) = Uo{x) = 1 and Ti{x) = x. Another useful recurrence relations are: 

(1-x^)T;;(x) = -nxT.a{x) + nTn^i{x) (A.2) 

Tyn{x)Tn{x) = ^ (Trn+n{x) + T\m-n\{x)^ m,n>0 (A. 3) 

The zeros of T„: Tn{C,j) = ; j = 1, 2, . . . , n are given by 

e, = cos^^:^^ ; j = l,2,...,n (A.4) 

The extrema points rjk at which |T„(a;)| = 1 are given by 

AjTt 

rjk = cos — ; k = 0,1,2, ... ,n (A. 5) 

n 

At these points we also have 

Tnivk) = i-lf ; k = 0,1,2,..., n (A.6) 

It is clear that the points rji, . . . , rjn-i are zeros of the derivatives T^{rjk) = ; k = 1, . . . ,n—l 
and therfore also zeros of Un-i{j]k) = ; k = 1, . . . ,n — 1. One also has 

T„(l) = l ; T„(-l) = (-ir (A.7) 

The Chebyshev polynomials take very simple form in the trigonometric representation: 

sinfra + l)p 

Tn\x) = cos np ; Un\x) = with x = cosp (A. 8) 

sin T 

One notices that for the special choice of arguments: 

P^ = ^^ ; j = 0,l,2,...,Ar (A.9) 
with some integer N the Chebyshev polynomials satisfy the periodicity relation 

Tn+N+lipj) = Tnipj) ; Un+N+lipj) = Unipj) (A.IO) 

Chebyshev polynomials become very useful in the study of the eigenvalue equation LijVj = 
Xvi of the matrix Lij related to the Cartan matrix. In components the generic equation 
reads as 

- Va-l + (2 - \)Va - Va+l = (A.ll) 

for the parameter a taking values between 2 and r — 1 = rank Q — 1. Introducing variable x = 
(2 — A)/2 into ([A.ll|) and rewriting linear combination of the Chebyshev polynomials 

Va= iATaix)+BUaix))vo (A. 12) 

we find that ( |A.11|) is trivially satisfied due to (|A.1|) . The remaining equations of KijVj = Xvi 
(whose not appearing in ( |A.11| )) impose the secular equation on the L matrix. They have 
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to be solved separately for each Lie group and solutions can be found in terms of the zeros 
and extrema of the Chebyshev polynomials as it was done in the text. 
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